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In this paper a set of approximate equations is derived which is applicable to very nonadiabatic, 
nondissipative, buoyant flows of a perfect gas. The flows are assumed to be generated by a heat source in which 
the heat is added slowly. The study is motivated by the occurrence of such flows in fires. There, the time scale 
associated with the fire growth and resultant fluid motion is usually long compared with the transit time of an 
acoustic signal (based on the temperature derived from the heat added) across the spatial extent of the fire. The 
appro- mate equations are characterized by a spatially uniform mean pressure appearing in both the energy 
equa on and the equation of state with the spatially nonuniform portion of the pressure only appearing in the 
momentum equation. Therefore, the pressure remains almost constant in space while significant density and 
temperature variations, such as might occur in a fire, are allowed. The approximate equations are shown to reduce 
to the Boussinesq equations when the heat addition is mild. These equations are also shown in general to admit 
internal-wave motions while "filtering out" high-frequency, acoustic waves. In addition, they are shown to be 
expressible in conservation form, the pressure satisfying an elliptic equation whose homogeneous terms are 
derivable from the wave equation by letting the sound speed become infinite. An equation for the mean pressure 
is also obtained. For the special case of a room heated at a uniform rate with a small leak to the outside, an 
approximate solution for the mean pressure is determined explicitly. 
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1 . Introduction 

The motions considered in this paper are those which arise solely due to localized addition of heat to an 
otherwise quiescent fluid in the presence of gravity. The principal application of interest to the authors is the 
movement of smoke and hot gases caused by fires. However, the formalism introduced will be appropriate to 
rederive equations previously used to describe rapid heating by a laser and equations which describe 
detonations, as well as those appropriate to fire research. 

The purpose of this paper is to obtain formally a set of equations of motion which permit description of large 
temperature and density variations due to volumetric heat addition without requiring a simultaneous 
description of acoustic oscillations arising because of the elastic properties of the fluid. Such model equations 
include the important features of buoyant flows without requiring excessive computer time necessary to 
determine high-frequency sound waves when numerically integrated. In this sense the equations "filter out" 
the sound waves while describing the lower frequency, organized motions due to buoyant effects such as 
internal waves. 

Such filtering is analogous to that employed previously to eliminate acoustic waves from the equations of 
motion in studies of atmospheric dynamics [1-4]. However, both the procedure used to implement the filtering 
and the final equations are quite different. Equations closely related to those derived in the present study 
have been used in attempts to describe both steady-state natural convection over a hot horizontal surface [5] 
and transient two-phase flow arising in nuclear reactor research [6, 7]. In each of these studies, large density 
variations due to temperature changes are admitted, but compressibility effects are suppressed. In the reactor- 
related papers, the liquid-vapor mixture is modeled as a thermally expanding, homogeneous fluid, and the 
discussions are closely tied to computer code development. In neither case has the effect of a closed 
environment on the pressure been considered in detail (see sec. 5). To the authors' knowledge, moreover, no 
systematic derivation of these equations has been presented previously. 

This derivation is different than that used when considering most atmospheric studies. Although the 
atmosphere is indeed thermally driven, the sources of heat are sufficiently diffuse that direct heating is 
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normally neglected, its principal effect being to stratify the atmosphere. The induced motions are in fact 
usually assumed to be nearly adiabatic disturbances to a pre-stratified state of rest (see ref . 1 , 8 and 9 for 
example). This stands in direct contrast to the nearly isobaric processes in any initially (nearly) unstratified 
environment of interest in fire research. 

Another difference is the role played by the Boussinesq equations. Considerable attention has been given 
in the geophysical and astrophysical literature to the conditions under which these equations may be used 
[10-12]. By contrast, their utility in fire research seems to the present authors to be rather limited, although 
many useful studies based on them have been performed. The equations derived below contain all the terms 
retained in the Boussinesq equations. However, the large variations in density caused by the local heat 
addition prevent the linearization of density about any nominal, unperturbed value. 

The fluid is taken to be an inviscid non-heat conducting perfect gas. The magnitude and the spatial and 
temporal variation of the heat source are taken as known. These approximations are justified because under 
conditions characterisitic of even a small room fire, the Grashof numbers (representing the ratio of the inertial 
to viscous forces for natural convection) are large enough for molecular transport phenomena to be important 
only in the highly convoluted flame sheets which constitute the region of intense heat addition. The study of 
the detailed flame structure of real fires is an extraordinarily complicated subject in its own right, and is 
bypassed here by specifying the heat source. Batchelor [16] gives a brief but relevant discussion of the 
applicability of the inviscid equations in the context of atmospheric motions. It should be noted that such 
simplifications do not preclude a description of turbulence; but no turbulence model is explicitly included in 
this study. Simple models of smoke and hot gas transport which neglect molecular transport phenomena have 
been reasonably successful in predicting global properties of flow fields [17]. The present work is intended as 
a first step towards more detailed studies along these lines. 

In the next section the model equations are derived using a dimensionless formulation based upon the 
characteristics of the (volumetric) heat source. In section 3 some of the properties of this system of equations 
are examined. Particularly, it is shown that internal waves will be described while acoustic waves are not 
permitted by this approximate set of equations. Also, the equations are rewritten in conservation form, a form 
often preferred for numerical computations, and the mixed hyperbolic, elliptic nature of the equations is 
noted. In the fourth section a discussion is given to show how the Boussinesq equations are obtained as a 
special case of these equations. In the final section, an equation for the mean pressure is obtained. This 
equation is novel and is a direct result of the filtering analysis performed in section 2. An approximate 
solution is obtained for the special case of an enclosure heated at a uniform rate with a small leak to the 
outside. 



2. Formulation 

The equations of motion for a perfect gas in the presence of an energy source of strength Q may be written 
in the form: 



— + — (pm) = 
dt dxi 

(du t duA dp 
PI — + "fc— + " Pg"i = (1) 

pC J?l + Uk dT\ _(§P + Uk §P 

\dt dxjcj \dt dxjc 

p = pRT 

Here p is the density, p the pressure, T the temperature, and u t the velocity at any position Xi and time t. 
Also, R is the gas constant, C p the constant-pressure specific heat, and n^g is the gravitational acceleration. 

The fundamental assumption on which this paper is based is that the source strength may be usefully 
characterized in the form: 
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Thus E is a dimensional constant which determines the strength of the source. The quantity Q is a function 
of order unity which varies smoothly with respect to its arguments. The independent variables X\ and t are now 
made dimensionless with respect to / and t respectively, where / and t are the length and time scales 
describing the spatial extent and the temporal variation of the heat source. The dependent variables are 
scaled as follows: 

P = p R(yu r) 



Po^v 

y- i£oV' 2 



p = 



y Po. 

y- 1 



Uyu r) 



(3) 
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The temperature scaling used in eq (3) is chosen so that the time derivatives in the energy equation always 
are the same magnitude as the source term. The density is normalized with respect to an ambient level p 
which occurs in the absence of any heat addition. The pressure scale then follows from the equation of state. 
The velocity is non-dimensionalized with respect to a thermal (sound) speed based on the temperature scale. 
The specific heat ratio is denoted by y. 

Substituting eqs (2) and (3) into eq (1); the dimensionless equations of motion are: 
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The parameter 8 is defined by 



5 = 



t (y - lEA* 



i\yp 0l 



This parameter plays a fundamental role in what follows. Physically, it represents the distance, based on the 

f y 



characteristic velocity 



1\ , \i 

Eq/ Po) -> a disturbance can travel during the time the source changes 



appreciably divided by the spatial extent of the source. The magnitude of 8 is directly related to the heating 
rate. The case 8 <^ 1 corresponds to rapid heat addition, realized in laser heating devices. The case 8 ~ 
(1), an intermediate case, describes detonation phenomena. Finally, the case 8 ^> 1, corresponding to slow 
heat addition, occurs in both controlled combustion and in fires. 

The mathematical task now reduces to displaying the dependence of the reduced variables defined in eq (3) 
upon 8 in the three cases described above. The dominant terms in the expansions of these variables as 
functions of 8 then satisfy equations which are subsets of eq (4). The expansions, and thus the resulting 
equations, are different for each case. They will now be treated in order. 
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(i) Fast heat addition (8 <^C 1). 
The dependent variables may be expanded in the form: 

P = P (0) ( ri ,r) + 0(8) 

R = R (0) (n .t) + 0(8) 

= #•> ( r , ,r) + 0(8) 

^=8{^ 0) (rj,T) + 0(8)}. 

Substituting the expansion (5) into eqs (4), the equations of motion reduce to: 



(5) 



dT 



= 



dr dy t E \y - 

'dd (0) y - 1 dP (0) ~ 



(6) 



#o) 



dr 7 dr 

p(0) — R(o)Q(o) 



Q(y t ,r) 



Equations (6) are a slight generalization of those used by Rehm [11] to study the laser bleaching wave. Note 
that once the heat addition model for Q is specified in terms of the thermodynamic variables and their 
gradients, the momentum equations are decoupled from the remaining three of eqs (6). Moreover, the density 
remains constant. This situation arises because the gas has not moved appreciably during the short interval in 
which the heat is added. The quantity p o gllE is usually quite small in laser applications. 

(ii) Intermediate heat addition (8 ~ 0(1)). 

When 8 is 0(1), no terms of eq (4) (except possibly those due to gravitation) can be ignored. In this case, it 



is easier to put 8=1. This is equivalent to replacing the length scale / by t 



y-iE Y 



Thi 



y Poj 

mathematical statement of the fact that the ratio of the length scale to the time scale is determined by the 
wave speed in a detonation [14]. The speed of such a front is of the same order as the thermal speed 

y-iE o y 
y po/ 

(iii) Slow heat addition (8 ^> 1). 

Since 8 is large, an expansion in negative powers of 8 is anticipated. Inspection of the energy equation 
shows that the spatial convection terms are too large by a factor 8 to be driven by the source term. If the 
velocity is reduced by a factor 8 _1 , the momentum balance requires the pressure gradients be reduced by a 
factor 8 -2 . The expansion then becomes: 



y. = g-i {y[o) + 0(5-2)} 

p = pio)(j) + 5-2p(l) 

= 0<°> + o(8~ 2 ) 

R = Rio) + ()(5-2) 



(7) 



The suppressed arguments of all dependent variables except P (o) (t) are both y t and r. Substitution of the 
forms (7) into eqs (4) yields: 



300 






l dr dy 5 J 
/? (0 M + ^ 0) \ 

I 5t a ri j 



^ — # (o) - o 

(8) 
y - 1 dP (o) ~ 
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R (o) e (o) = p(o)( T ) 

Note that the expansion (7) implies that the spatial variations in pressure are extremely small. The 
reference pressure level ^^(r) can be determined from the boundary conditions. For example, the pressure 
variation in an enclosed volume W subject to heating is given by (see sect. 5 for a more detailed derivation): 

1 dl* 0) f ~ 

-v-r=\ Q d3 r- ( 9 ) 

y dr J w 

If the domain is open to the atmosphere, a point can be found at which the pressure is specified as Pa,, say. 
Then: 



fXo 



1 E n 



This reference pressure point can be taken as the origin of the vertical coordinate y^ without loss of 
generality. 

If the quantity gt%/l <^ 1, then the gravitational force may be neglected. This situation arises frequently in 
studies of nearly isobaric combustion processes [13]. However, in fire research there are many problems for 
which a = gt 2 /l is 0(1). It is then convenient to redefine the time scale t such that a = 1. The velocities are 
then 0((g/)*) in magnitude, while the spatial pressure perturbation is 0(pgl) and the temperature is 



.fi-l 
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In dimensional form, the approximate equations of motion (eqs (8)) are: 

d P * • ^ 

— + — (p Uj ) = 

dt dxj 



(dui duf 

dt dxj 



(11) 



, d(p ~ pa) 
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( dT df\ d Po 

pC P [ — + u 5 — - — = Q 

H \ dt 3 dxj dt V 

Po (t) = pRT. 

Since p is necessarily independent of position, it follows from the equation of state that the ambient 
atmosphere in the absence of heating must be taken as having a constant density. This limitation is ultimately 
traceable to the assumption stated earlier that the vertical extent of the source is much less than the scale 
height of the atmosphere. It is often convenient to replace the variable p — p by a reduced pressure /?*, 
defined by: 

p* = p - Po~ PogXiKi- (12) 
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Equations (11) are then altered in that p must be replaced by p — p in the gravitational force term and 
p — p replaced by/?* in the pressure gradient term. With this transformation, eqs (11) appear more similar 
to the Boussinesq equations. 



3. Properties of the Derived Equations 



In this section two properties of eqs (11) will be examined. First, it will be demonstrated that these 
equations contain buoyant effects without including the elastic effects of the medium. This is demonstrated by 
examining small-amplitude waves and showing that internal waves (buoyant effects) are allowed in this 
description while acoustic waves (elastic effects) are excluded. Second, the conservation form of eqs (11) are 
examined and the mixed hyperbolic and elliptic nature of these equations is noted explicitly. 

To demonstrate the first property, we assume that the flow has evolved to point where, in eqs (11), the 
mean thermodynamic variables are stratified so that, p = p 00 {z) and T = T 00 (z) while/) = p = constant. (For 
convenience, the spatial coordinates will be denoted by x, y and z.) All the flow variables are perturbed about 
these quantities. For simplicity, the heat source is taken to be zero. Also, we will take the buoyancy frequency 

g J p 00 to be constant (so that the density is stratified exponentially). The type of waves described 



N l 



2 -- _ 



dz 



by eqs (11) is not changed by these simplifications. The equations (with variable N 2 {z)) would describe the 
fluctuations produced in a closed room following a fire which had stably stratified the environment outside the 
buoyant plume and then gone out after consuming the available oxygen supply. 

The perturbed dependent variables are substituted into eqs (11), and the equations are linearized in the 
perturbation quantities. The temperature perturbation can be eliminated from the resulting equations. 

From the momentum equations and the continuity equation, the velocities can then be eliminated to give an 
equation coupling the density and pressure perturbations. 

From the energy equation and ^-component of the momentum equation, the vertical velocity can be 
eliminated to give a second equation coupling the density and pressure perturbations. These two coupled 
equations can then be used to obtain a single equation for the perturbation pressure p: 



d 2 d 2\ d 2 d 2 tf d 2 d 
+ N 2 + + + 

/ \ dx 2 dy 2 / dt 2 dz 2 g dt 2 dz 



The dispersion relation is obtained by assuming waves of the form 



p=0. 



[- 
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p = A exp z + i(xk x + yk v + zk z — cat) 
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(13) 
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Then 




or 



(15) 



This dispersion relation is the one which describes internal gravity waves. For the two-dimensional case 
(i.e. when there is no dependence of the motion in the y-direction), the equation was derived and the waves 
discussed, for example, by Lamb [18], by Mowbray and Rarity [19] and by Whitham [21]. The equations do 
not permit acoustic waves, as claimed earlier. In the derivations referenced above, the starting equations are 
taken to be those appropriate for an incompressible fluid. The dispersion relations and properties of more 
general acoustic-gravity waves in a compressible, stratified medium are discussed in books by Eckart [8], and 
by Yih [20] and by Dutton [9] for example. 

It is important to note that the model equations contain only internal waves and not acoustic waves. 
Apparently some confusion has existed in the past concerning what assumptions are necessary to assure this 
filtering. For example, in the paper by Gough [4], where the "anelastic approximation" is derived and 
discussed, the continuity equation is replaced by V*(pu) = 0, and it is stated that this approximation is 
necessary in order to remove the acoustic modes of oscillation. As seen above, the time derivative of the 



302 



density in the continuity equation need not be eliminated, however. Rather, acoustic modes are removed 
when the pressure in the momentum equation is decoupled from density and temperature fluctuations arising 
through the equation of state. 

Many numerical methods for solving fluid mechanics problems are based on finite difference schemes 
obtained from conservation forms of the equations of motion. The appropriate forms of the continuity and 
momentum equations are given by: 



— (pud 
at 



d P d , X 

— + — (pu t ) = 
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(16) 
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Equation (16) may be thought of as allowing the computation of p and U{ at a time t + 8t given information 
about all quantities at t. In order to complete the hypothetical computation of the flow variables at t + 8/, 
equations for the pressure must be determined. This determination consists of two parts; an equation for the 
reference level p (j), and an equation for the spatial distribution of the overpressure p(x h i). The equation for 
p {t) is derived and discussed in section 5. A solution for this reference level pressure is also presented. 

The equation for the overpressure/? is obtained by dividing the momentum eq (11) by p and taking the 
divergence of the resulting expression. Then, using the continuity equation: 
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dxi 

Note that the terms not involving/? in eq (17) are all known at time t + dt. Since the density p is also known 
at t + 6t, p satisfies a linear self adjoint elliptic equation at each instant of time. The appropriate boundary 
conditions are a specification of the normal derivative of/> obtained from evaluation of the momentum equation 
at the boundary. Such a specification ensures that the consistency condition for the existence of a solution to 
eq (17) is satisfied automatically. This equation is the generalization of Poisson's equation appropriate to an 
inhomogeneous fluid. It remains an elliptic equation, despite the density variations, because the sound waves 
which would lead to a finite propagation speed (and hence to a hyperbolic system of equations) have been 
eliminated. 

This can be seen quite readily if the steps leading to eq (17) are performed on the full equations (eqs ( 1 )) 
rather than the reduced set (eqs (11)). The result is: 
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dt 
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The quantity C 2 in eq (18) is the square of the local sound speed. Equation (18) is the wave equation for the 

pressure, a hyperbolic equation generalizing eq (17). The terms in the second brackets on the right hand side 

1 d 2 p 
above, together with the term — — — - , also appear on the right hand side of eq (17). 

Equation (18) can be reduced formally to eq (17) by nondimensionalization and expansion of the dependent 
variables (retaining only terms 0(1) and neglecting 0(8 -1 )) as in section 2. When this is done, terms linear in 
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pressure gradients and quadratic in velocities must be retained. Products of pressure gradients and velocities 
are of higher order and may be ignored. Then eq (18) reduces to eq (17). 

The overall system of equations is neither elliptic nor hyperbolic. The mixed character of the equations 
may be seen by noting that the energy equation (eq (11)) is a hyperbolic equation for T. Indeed, the overall 
system may be either elliptic or hyperbolic (or neither) depending on the problem under consideration. 

The variable density coefficient in the pressure equation complicates the task of numerical solution. 
However, relatively efficient schemes for handling equations of this type have been developed [22]. Thus, 
practical finite difference computational schemes based on eqs (16), (17) and an equation for the mean 
pressure (discussed in sect. 5) are feasible. 

4. The Boussinesq Equations 

The derivation presented in section 2 is based on two assumptions. First, the length scale, time scale and 
the temperature scale associated with the volumetric heat source are such that the heat addition is slow. This 
assumption implies that the pressure over a large region surrounding the source is almost uniform in space 
(while it may vary with time) during heating. However, it does not imply any restriction upon the magnitude 
of the density (or temperature) variation during heating. The second assumption is that the flow velocities are 
induced by buoyant effects. This assumption relates the magnitude of the temperature variation, the density 
variation and the flow velocities induced by the heat source. 

Batchelor, in obtaining conditions for which atmospheric motions are determined only by an overall 
Richardson number (the ratio of buoyancy to inertia forces), has discussed the possibility of atmospheric 
motions where the pressure varies only slightly from hydrostatic, while the density and temperature deviate 
significantly from hydrostatic values [16]. Batchelor was concerned explicitly with adiabatic flows. However, 
for nonadiabatic motions the same possibility exists when heat is released slowly, as shown by the derivation 
in section 2. 

The Boussinesq equations arise when two additional assumptions are made. The first assumption is that the 
scale height associated with the static density variation (the density variation in the absence of motions) is 
much larger than the vertical length scales of interest. Then the static density variation from its mean value is 
small. The second assumption is that the density variations produced by the heat source are small. These are 
the assumptions made by Speigel and Veronis [10] in their derivation of the Boussinesq equations. 

The height scale associated with hydrostatic variations in the thermodynamic variables was already 
assumed to be much larger than the vertical length of interest in section 2. Therefore, the spatial variations in 
temperature and density (as well as in pressure) in the absence of motion were taken as small, satisfying the 
first assumption stated by Speigel and Veronis. (In cases of practical interest in fire research, the hydrostatic 
variations of these quantities are negligible compared with variations induced by the heat source.) The second 
assumption, that motion-related temperature and density variations are small, implies that the temperature- 
variation scale is small compared with the nominal ambient temperature T : hence we assume = 

Po^p Po^v' o 

With these additional assumptions (and a change in the time scale discussed below), the Boussinesq 
equations can be derived formally from eqs (8) by expanding the dimensionless dependent variables in a 
series of powers of e. (This second expansion implies that the Boussinesq equations are obtained by a two- 
parameter expansion of dependent variables in terms of parameters 1/8 and e.) However, rather than 
performing this second expansion, which is very formal and notationally cumbersome, we demonstrate how 
the Boussinesq equations are obtained from the dimensional equations, (11), in a more physically meaningful 
way. 

As suggested in section 2, we introduce the reduced pressure p* defined in eq (12). We also define a 
reduced density p* and a reduced temperature 7 1 * as follows: 

P* = P ~ Po 

(19) 
V = T - T . 

For simplicity we assume thatp is constant. Introduction of definitions (12) and (19) into Eqs (11) yields 
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(20) 



Po = (p. + p*)(r + 7*)«. 

The assumption that motion-related temperature and density variations are small implies that the reduced 
density, temperature and pressure must be proportional to e. Since the buoyant effects drive the flow, the 
acceleration term in the second of eqs (20) must also be proportional to e. This implies that the time scale for 

motion is t' — tje^ = I — I , and that the velocity scale is proportional to 6*. Finally, since Q « E /t , the 

w 

heat source is proportional to €*, Therefore, with the formal substitutions 

P* = ep 

p* = ep 



— ^2// 



T* = eT 

U = €^U 

t = t'/ei 



Equations (20) become, to lowest order in €, 



du* 

— = o 

dxj 

(dui A duA dp 
at oXj/ oxi 

(22) 
(of A dt\ 

P f 

— + — = 
Po T 

where p = RpJ . 

These are the equations, in the absence of dissipative effects, usually described as the Boussinesq equations: 

the density is considered to be constant except where the density difference produces a buoyant force [10]. 

5. The Mean Pressure 

The variable mean pressure level p (t) is determined from the approximate system of eqs (11) by a 
consistency argument. The continuity equation (the first of eqs (11)) may be multiplied by C P T and added to 
the energy equation (the third of eqs (11)) to obtain: 

PoW T~ + --— = Q- (23) 



y — 1 dx k y — 1 dt 
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As mentioned earlier if the domain under consideration is open to the ambient atmosphere, thenp is constant 
and need not be considered further. Now let the domain be a "closed" volume V (the meaning of "closed" will 
be made clear) bounded by a surface A. Integration of eq (23) over V yields: 



— - I UiUi i 



dA \ Po (t) + 



V d Po 

y - 1 dt 



-I 



QdV. 



(24) 



Equation (24) is a differential equation for p (t) whose coefficients are determined by the geometry of the 
enclosure and boundary conditions specifying the normal component of the velocity on the boundary. Thus, 
p (t) is determined by the requirement that the net efflux of fluid from the enclosure be consistent with the net 
rate of heat addition. Note that if the room were truly closed, in the sense that ufii = at each point on the 
boundary, the mean pressure would rise linearly for a constant heating rate. This would very quickly lead to 
pressures sufficiently high to cause structural collapse. Most rooms in buildings are not closed in this sense, 
even with doors and windows shut. A more realistic geometry will be discussed below. 

Consider a room of volume V whose general shape is like that illustrated in the insert in figure 1 . The room 
is open to the ambient atmosphere at pressure p^ only through a crack of length / and height d. The aspect 
ratio lid and the Reynolds number of the crack are assumed to be sufficiently large for the flow in the crack to 
be effectively two dimensional and inviscid. The crack is taken to be near the floor but away from the fire. 
The density of the gas exiting through the opening will then be the ambient value p^ until the hot ceiling layer 
formed by the fire plume descends to the floor. Under these circumstances, the efflux is given by the usual 
orifice formula [23] 



UiUi dA = — — ld\ — (p {t) - poo) f • 
J 77 + 2 [poo J 



(25) 



FORCE ON 2 METER SQ. 
OFFICE WINDOW 
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FIGURE 1. The force on a two square meter window caused by the mean pressure rise in a room during afire. The fire sizes correspond 
to a small and a large waste-container fire while the room dimensions are those of a typical N.B.S. office. The room is assumed to be closed 
except for a crack of dimensions 75cm by 2.5cm. 



The total rate of heat release for various burning materials is a frequently measured parameter in fire 
research. A useful compilation is given in ref. 17 and the references therein. For the present calculation this 
quantity will be characterized as follows: 



/ 



QdV= Q 



(26) 
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The quantity Q will be taken as the total rate of heat release due to the fire. The losses due to heat transfer to 
the boundaries are unimportant on the time scale associated with the pressure buildup, although ultimately 
they are a substantial fraction of Q. 

It is convenient to introduce non dimensional variables forp and t. 



PoU) = Ml + e 2 P(r)} 

m y - 1 (ir + 2\ Q 



V2y \ ir J I dCnP 
tt + 2\ 2 (y - l)QV 



(27) 
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The quantity Coo in eq (27) is the sound speed based on the ambient pressure and density. The parameter € 
determines the maximum steady state pressure rise attainable for a given fire size and enclosure geometry. 
For most realistic scenarios, e < 1. Substitution of eqs (25)-(27) into eq (24) and ignoring terms 0(e 2 ) leads to 
the evolution equation for p(r) 



dP 
dr 



1 - y/P. (28) 



The solution is readily obtained in the form 

T = log (1 - V*V - y/P. (29) 

Note that P = 1 corresponds to the steady state, so that the final fractional pressure rise measured in 
atmospheres is equal to e 2 . Thus, the pressure rise is directly proportional to Q 2 and inversely proportional to 
the square of the opening area. Although eq (29) can be plotted as a single curve, it is more instructive to pick 
a specific geometry and display the results in physical units as in figure 1. The fire sizes chosen correspond 
to a small and a large waste container fire [17], while the room dimensions are those of a typical NBS office. 
The crack dimensions are taken to be 75 cm by 2.5 cm. The quantity e 2 is approximately 0.01 for the larger 
fire and the Reynolds number based on crack height is 0(10 4 ) for the smaller fire. 

At a time t = 2.5 s only a few percent of the original room air has left through the crack, so the hot ceiling 
layer is unlikely to be very deep. Thus, the assumptions underlying the calculation are amply fulfilled. It 
seems that even a rather unspectacular fire can rapidly generate a considerable unbalanced force on a 
structure, which can rupture a relatively weak component like a window. This may be at least a partial 
explanation of the fact, commonly observed, that windows often break in fires. 

In summary, the equations derived in section 2 are ones which we feel are important in fire research and 
may be important in other applications as well. The flows considered are highly nonadiabatic and buoyancy 
effects are dominant, yet the pressure remains almost uniform spatially as the flow evolves. This implies that 
the time scale associated with the heat source, which is assumed to be comparable with the time scale for 
buoyancy effects, is much larger than that time scale required for acoustic signals to equilibrate the pressure 
over the region of interest. The former time is the one of interest for such flows, and therefore, it is desirable 
to remove the acoustic waves so that the time step size for numerical integration of the appropriate equations 
is not limited by the time scale associated with these waves. The equations derived in section 2 accomplish 
this filtering, as shown in section 3. Furthermore, the equations can be written in conservation form, a form 
often preferred for numerical computations. 
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